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Abstract 

We propose an alternative way of looking at data on anomalous 
J/ip suppression. The proposed method is in principle equivalent to 
the one used by the NA50 Collaboration, but it permits to visual- 
ize separate contributions of individual processes responsible for the 
disintegration of J/^'s produced by a hard process in nuclear colli- 
sions. The method can be used provided that the time sequence of 
contributing mechanisms is known or assumed. It offers an alternative 
graphical presentation of the onset of anomalous J/ip suppression in 
Pb-Pb interactions observed by the NA50 Collaboration at the CERN 
SPS and might contribute to explain why different mechanisms, such 
as J/ij) suppression by the Quark-Gluon Plasma and by co-movers in 
the Dual Parton Model or in Monte Carlo microscopic approaches, all 
lead to an approximate description of anomalous J/ip suppression. 



1 Introduction 



The anomalous J / ip suppression observed by the NA50 Collaboration |], || in 
Pb-Pb collisions at the CERN SPS can be naturally interpreted as J/ip disso- 
lution by the Quark-Gluon Plasma (QGP). This phenomenon was predicted 
by Matsui and Satz |3j more than 15 years ago. The interpretation of the 
anomalous J /ip suppression is somewhat complicated by the fact that there 
are also other contributions to J/ip suppression. Nuclear absorption (NA) 
(sometimes referred to as "pre-resonance absorption" or " Gerschel-Hiifner 
mechanism") 0-0 is probably responsible for most of Jjip suppresion in 
nuclear collisions induced by lighter ions. Disintegration of J / ip by collisions 
with secondary hadrons ||7||— |[L3| (also referred to as "co-mover" interaction) 
is another possible contribution. In order to make the J/ip suppression by 
the gas of secondary hadrons sufficiently large, in some approaches [fL0| , [12] 
cross-section (<t co ) for the process hadron-l- J/ip — > DD + X has to be of 
about 5 mb, depending also on the value of the cross-section (a a ) for nuclear 
absorption. In the approach based on the Dual Parton Model (DPM) [|TTJ 
one needs a co = 0.6 mb, for a a = 6.7 mb and |T3j a co = 1 mb for a a = 4.5 mb. 
It has been argued [14] that, at energies corresponding to the thermal motion 
of hadrons, the cross-sections for the disintegration of J / ip by hadrons are in 
fact an order of magnitude smaller. The issue is not yet definitely clarified. 

There is one key difference between the suppression of J/ip by the QGP 
and by other mechanisms. For J/ip dissolution by QGP a rapid onset of 
the anomalous suppression is expected or at least possible, whereas for other 
mechanisms one expects a smooth dependence of J/ip suppression on the im- 
pact parameter b and on the total transverse energy of a nuclear collision. 
Models by Blaizot and Ollitrault (BO) |TH] and by Kharzeev, Lourengo, Nardi 
and Satz (KLNS) [16|, [HJ do show a rather rapid onset of the dependence 
of the anomalous J/ip suppression on E? or on b. Plots of experimental 
data on J/ip suppression, in particular those with the nuclear absorption 
"subtracted" |T3| do show such an abrupt onset. 

Although the anomalous J/ip suppression is a rather spectacular Phe- 
nomenon, indicating a rapid onset of a new mechanism of J / ip Suppression, 
there are a few mechanisms that are able to obtain some agreement with the 



basic features of the data. Apart from the suppression by QGP [1151 16, 17 



these include the disintegration of J / ip by co- movers JTl], [L3| and Monte Carlo 



microscopic models [10|. In order to see why it is possible that rather differ- 
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ent models lead to a similar overall J/ip suppression, it would be interesting 
to see not only the resulting total J /if) suppression but also contributions 
of individual mechanisms involved in this total suppression. The purpose of 
the present paper is to propose an alternative way of looking at the data on 
J /if) suppression, and in particular on the onset of anomalous J /if) suppres- 
sion. Our approach is based on the additive decomposition of contributions 
of different mechanisms responsible for the disintegration of J/ip's produced 
by a hard process in nuclear collisions. The method can be used provided 
that the time sequence of contributing mechanisms is known or assumed. In 
some cases there are good reasons to believe that two mechanisms are at 
work in the same time period. This would concern, for instance, the nuclear 
absorption and possible depletion of gluon structure functions |TI| during the 
first part of the nuclear collision. In such both mechanisms should be 

considered as one stage. 

In the next section we shall describe the method; in Sect. 3 we shall 
present a few illustrative examples. In Sect. 4 we shall discuss the relation- 
ship of the proposed way of looking at the data with the standard procedure 
used by the NA50 Collaboration, and comments and conclusions will be pre- 
sented in Sect. 5. 



2 Separation of different contributions to J/ip 
disintegration 

We shall consider here for simplicity the case when only two mechanisms of 
J /if) suppression are present: the nuclear absorption and the dissolution of 
J/ip by QGP, the latter taken from a point of view very close to that of the 



BO H and KLNS 0, models. 
We shall use the following notation 

Nplo d (E T ) = the total number of J /if) produced in collisions at given Et-i 

Nha(Et) = the number of J/ip's disintegrated by NA, 

Nqqp(Et) = the number of J/ip's dissolved by QGP, 

n S( e t) = the number of surviving J/ip's, measured in experiment 

N^ er (E T ) = the numbers of J/ip's disintegrated by other mechanisms. 
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The balance between input and output requires 



nZ( E t) = Ng(Br) + N J Jt{E T ) + N J ^ P {E T ) + N^ er (E T ). (1) 



We suggest that all terms in Eq. (1) be calculated and plotted separately, 
or in combinations, for a given A-B collision, for all values of Et- In this 
sense our approach tries to give an answer to the question "Where are all the 
J/^'s gone?" or corresponds to "an accountant's look at the data on J/ip 
suppression" . 

Admittedly, the suggestion we are making is rather trivial, but we do hope 
that looking at the data in this way can make the results of a particular model 
more transparent and can help to avoid some inconsistencies. 

In specific applications Eq.(l) can be rewritten in such a way that more 
reliably known terms are put on one side, thus providing constraints for the 
more model-dependent and less reliably known ones. We can for instance put 
on one side the best known terms N^ d (E T ) — (E T ) or even N^ d (E T ) — 

N eip( E T) ~ N J Jf{E T ), and try to describe or fit these expressions by a 

particular model. For instance in the latter case by Nqq P (E t ) + N^ er {E T ). 

The purpose of the present paper is not to attempt a detailed analysis of 
the data. We wish just to describe the method and to give a few examples 
as illustrations. For that purpose we shall take nuclei as homogeneous hard 
spheres with radii Ra = 1.2A 1//3 fm, we shall neglect the energy fluctuations, 
including the "knee" in multiplicity distribution at E T xs lOOGeV and we 
shall assume here a strong correlation between b and E T : 



where N w (b) is the number of interacting ("wounded") nucleons. With these 
simplifications the quantities entering Eq. (1) can be expressed in the fol- 
lowing way: 



nucleon collision and N co u(Et) is the number of nucleon-nucleon collisions 
at a given value of Et- In a rather simplified case we have 



E T (b) = 0.325GeV N w (b), 



(2) 
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f R A sds f 2n c l a(s) rL B (b,s,e) 

/ / dO / PA°nndz A \ PB&nndZB, (4) 

Jo <T nn Jo J-L A (s) J~L B (s,e) 



where 



L A (s) = yjR\ - s 2 , L B (b, s, 9) = yj R 2 B — b 2 — s 2 + 2bscos(9) , 

when the expression under the square-root in Lb is negative, there is no 
tube-on-tube collision and the contribution vanishes. Nuclear densities are 
denoted as Pa, Pb (we are working in the approximation of nuclei as hard 
spheres, therefore Pa,Pb are constants), b is the impact parameter, s is the 
distance from the centre of the A-nucleus and 9 is the angle between b and s. 
A nuclear collision is taken as the sum of tube-on-tube collisions, with lengths 
of tubes 2La and 2Lb, and z a ,zb specifying the coordinate of nucleons within 
both colliding tubes: za varies from — La to La, zb from — Lb to Lb, and 
both za and zb increase in the direction of motion of the tubes in the cm. 
frame of the nucleon-nucleon collision. The nucleon-nucleon non-diffractive 
cross-section a nn is taken as 30mb. Experimental results on the survival 
probability S(E T ) of J/ip enter the expression 

NtH{E T ) = S{E T )N J p 't d {E T ). (5) 
Nuclear absorption of J/ip is given by the expression 

jii [Ra sds f 2n C^a{ s ) rLB(b,s,9) 

N NA ( E t) = elln \ / d0 / PA^nndz A / pB^nndz B 

JO a n n JO J~L A (s) J-L B (s,e) 

_ e -PAVa{zA+LA(s)} e -pBVa{z B +L B (b,S,0)}^ ^ 

where a a is the cross-section describing the absorption of J/ip by nucleons. 

The integral in Eq. (6) can be decomposed into two separate integrals, 
the former giving the total number of J/^'s produced and the latter being 
equal to the number of J/^'s that survive the nuclear absorption. The term 
Nqqp(Et) gives the number of J/^'s destroyed by the QGP; it is given as 



ji i fRA sd-S f 2n cLa(s) fL B (b,s,8) 
N QGp( E t) = Onn \ / d9 pA<?nndz A / pB&nndz B 

Jo a nn Jo J-l a (s) J-L B ( s ,e) 



r R A Sds f 27r ,„ /" L a(s) fL B (b,s,8) 

PA^nndZA I 
-L A (s) J-L B (s,d) 

e -PA<Ta[z A +L A (s)} e -p B a a [z B +L B (b, S ,0)]Q( K _ KcHt y (7) 
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Following KLNS [|IJ, see also Ref. 



, we have introduced 



k _ PA^nn^L A (s).p B a nn 2L B {b 1 s, 9) 
p A a nn 2L A (s) + p B a nn 2L B (b, s, 9) ' 

The parameter K crit specifies the onset of the QGP formation. This parameter 
has to be determined by the data in such a way that 

N J p H d {E T ) = N%f{Er) + NiH{E T ) + N J Q ^ P (E T ), (9) 

where we have assumed that the contribution N^ er can be neglected. Equa- 
tion (9) simply says that out of all J/^'s produced some survived, others 
were disintegrated by the nuclear absorption and others by QGP. 

The multiplicative factor is present in all of the terms in Eq. (9) 
and in what follows we shall leave it out. 

The time ordering of the different processes enters Eq. (7). It is most 
likely that the QGP is formed after the colliding nuclei have passed through 
one another and the QGP can thus destroy only those J/^'s that survived 
the nuclear absorption. 



3 Illustrative examples 

Figures showing the experimental data in the way they were used by NA50 
Collaboration [| are well known. In these figures they plot the survival 
probability of J/ip as 

0/771 \ J /ip measured / 1 n\ 

S{Et) = j/vw^, ' (10) 

inferred from direct data on J/ip over Drell-Yan pair production. The co- 
efficient S(Et) contains in a multiplicative way probabilities that J/ip has 
survived nuclear absorption, QGP dissolution and other possible mechanisms 
of J / ip disintegration. In the way of looking at data proposed here we plot 
separately the surviving J/ijj's (as given by experiment), those disintegrated 
by nuclear absorption, and those dissolved by the QGP and possibly also by 
other mechanisms. 

In Fig. [l] we plot as an illustration the terms N^ d , (survivors), and 

Nft A: leaving out the common factor o J J^ in Eqs. (3)-(7). The anomaly is 
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Figure 1: Results on N'^ d as calculated by Eq.(3) with a nn =3fm 2 , the data 
of NA50 Collaboration g, | used to calculate Ngg by Eq. (5) and the N J Jf 



calculated by Eq.(6) with a a =0.7fm 2 , and Nqq P evaluated by Eq. (7) 



not easily visible, since it appears as a rapidly growing difference N^ d —N^. 
In order to make the anomaly visible, we plot in Fig. |2] the expression AiV 
= ^prod~ ^eip ~ ^na with a nuclear absorption cross-section a a = 0.7 fm 2 . 
The shape of AiV indicates a presence of a J/if> dissolving mechanism with 
threshold between E T = 30GeV and E T = 40GeV. Note that values of AiV 
are slightly negative for the lowest three E? points, which is connected with 
the fact that in the standard presentation of data calculations based on the 
NA mechanism are below the lowest E? data. 

In Fig. |3| we show the opening of the space for other contributions, most 
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Figure 2: The expression AN = N^J d - N^ff - N J Jl as a function of E T . 
The term iV^f has been calculated with a a =0.7fm 2 . 

probably for J j ip suppression by a hadron gas, when a a in nuclear absorption 
becomes smaller. 

The true consistency check of the phenomenological description of data 
is provided by Eq. (1) with given by Eq. (5) and S(Et) taken from 

experimental data on J/ijj survival. In the illustrative case discussed above 
we have N'^ er = and Eq. (1) reduces to 

N J p H d {E T ) = E(E T ) ee N^(E T ) + N J Jt{E T ) + N^ P (E T ). (11) 

In Fig. |] we compare N^^Et) and H(Et)- The agreement is quite 
reasonable, in view of simplifying assumptions that we made. Agreement 
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Figure 3: The expression AiV = N^.^ — — iV^f as a function of 

Et- The term has been calculated with o a =0.7fm 2 , a a =0.6fm 2 , and 

<t„ =0.5fm 2 . 



with data is worse for the points with the highest values of Et where a 
large part of Ex is given by fluctuations |T3[ and our simplified model is 



not applicable. Note that the agreement visible in Fig. £| is non-trivial since 
N^p(Et) = S(ET)Nplg d is given by experimental data on S(Et). 

The information contained in Fig. [I] and Fig. ^ can be compressed into a 
single Fig. §, where we plot N J Jt(E T )- N J Jt(E T )+ N^ P (E T ) and N^(E T )+ 

Nqq P (E t )+ N^(E t ); and for comparison with the last expression also 
N J A{Et). 
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Figure 4: The comparison of iV ^ (solid line) and E(-Et) (open circles). 

4 Common points of the new way of looking 
at data and the standard procedure 

We shall first show that the way of plotting experimental data and results of 
phenomenological calculations is in some aspects equivalent to the standard 
way introduced and used by the NA50 Collaboration. 

To make the argument transparent, suppose that a certain number Nq 
of J/V's is passing successively through three obstacles (media) and can be 
disintegrated by each of them. The probability that it passes (survives) 
through the first obstacle is denoted as P\. The number of J/^'s passing 
through is N Pi, the number of those that fail to pass is Npi = N (l — Pi). 
In the same way the number of J/^'s passing through the second obstacle 
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Figure 5: The compressed information from Figs. [I] 

and | N J Jt(E T ) (dash-dotted line); N J J%{E T ) + N J ^ P (E T ) (dashed line); 

N J Jt{E T )+N J '£ P {E T )+N J e /${E T ) (open circles) and N^ d (E T ) (solid line). 



is N0P1P2 and the number of those that passed the first obstacle, but failed 
to pass the second one is N F2 = NqPx(1 — P2). Finally, the number of J/ip's 
passing throught the third obstacle is N PiP 2 Ps and the number of those 
that failed to pass the third obstacle is N F3 = N P 1 P 2 (l — P 3 ). The obvious 
identity 

N = N (l - Px) + N Pi(l - P2) + N P 1 P 2 (1 - P 3 ) + PxP2P.N0 (12) 
can be rewritten as 

iVo = N F1 + N F2 + N F3 + SN . (13) 
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where 

S = Pl P 2 P 3 = N ° ~ N ™ ~ N * ~ ^ 3 (14) 

iVo 

gives the overall survival probability. When comparing the results of the 
calculations with the data in the usual way the comparison is between exper- 
imentally measured values of S and the probability PiP 2 P 3 calculated from 
phenomenological models. The alternative way used here consists in proceed- 
ing according to Eq. (12), plotting individual terms N F1 , N F2 , N F3 , SN , with 
S taken from experiment and verifying the validity of Eq.(12). 

In a real situation, each J/ if) is born at a certain value of the parameters 
s, 9 (when considering nuclear collisions at a fixed value of b, and at a cer- 
tain value of b, s, 9 when integrating also over b. The probability of passing 
though different media depends on values of these parameters. Equations 
(11-13) change only a little. Lumping s, 9, or b, s, 9 into a single parameter x 
and introducing probabilities Pi(x), P2(x) and Ps(x), we have in an obvious 
notation: 

N = [n (x)dx 



N P1 = J n (x)P 1 (x)dx, N F1 = J n (x)[l - P x (x)\dx 
N P2 = J n (x)P 1 (x)P 2 (x)dx, N F2 = J n (x)P 1 (x)[l - P 2 (x)]dx 

N P3 = J n {x)P 1 {x)P2{x)P 3 {x)dx, N F3 = J n (x)P 1 (x)P 2 (x)[l - P 3 (x)]dx 

and one can write again Eqs. (11)— (13). 

When calculating the total J / ip suppression via an expression correspond- 
ing to Eq. (14), we get: 

M J/if> M J/i> A tJ/4> 

S ( E T) = -J7^ • (15) 

prod 

Using Eqs. (4)-(6) and the identity 1 — Q(k — K crit ) = Q(n cr u — k), we obtain 
from Eq. (15) the standard expression valid for J/ip suppression by nuclear 
absorption and by the QGP 

. fRA sds f 27T [ L a{s) rL B (b,s,6) 

S '' V (E T ) = / d9 pA<?nndz A / pB<?nndz B 

JO a nn Jo J-L A (s) J-L B (s,9) 

e -PA«a{z A +L A (s)] e -p B a a {z B+ L B (b,s,e)]Q^ crit _ K y ^ 
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5 Comments and conclusions 



We have desribed here an alternative way of plotting the data and model 
calculations on J /if) suppression and presented a few illustrative examples. 
A more detailed analysis of J / if) suppression along this path would certainly 
require using more realistic Woods-Saxon nuclear densities, a more realistic 
relationship between b and E T , detailed analysis of J /if) suppression in colli- 
sions induced by lighter ions at 200GeV per nucleon and a detailed evaluation 
of experimental errors. We are of the opinion that such an undertaking, al- 
though rather tedious, might help us to understand why models based on 
rather different assumptions lead to roughly similar results. 

For models considering only nuclear absorption and J / if) interaction with 
co-movers, as in Refs. [ll], [12], 13] the present scheme requires only minimal 



modifications. Equation (9) should be rewritten as 

KltAEr) = N J & ${E T ) + N J Jt(E T ) + N^(E T ). (17) 

where the last term corresponds to the number of J/^'s suppressed by the 
interaction with co- movers (with the gas of secondary hadrons). The term 
N^ d (E T ) is given by Eqs. (3) and (4), N^(E T ) by Eq.(5), N J Jt{E T ) by 
Eq. (6) and N^J^(Et) is calculated, in analogy to Eq. (7), as 



T/ . T/l [Ra sd,S f 27T C L a(s) rL B (b,s,6) 

N^(E T ) = <T J J* / dO PACTnndZA / p B (T m dz B 

J0 Onn JO J-La(s) J-L b (s,9) 

e -PAa a [z A +L A (s)] e - PB a a [z B +L B (b,sM _ gxp (_ f dr(va co p co ))) . (18) 

In Eq.(18), a co is the cross-section for the disintegration of J/if> by the inter- 
action with co- movers, p co (t) is the density of co- movers as a function of the 
proper time r of J /if). For more details see Refs. ||, |], [TTj, [12], [13, |2T |. 



an an- 



For Monte Carlo microscopic models such as the one in Ref. JTTJ 
alytic recipe for separating different contributions cannot be written down, 
but the separation is most likely possible as well. 

The alternative way of plotting and looking at data might be useful, 
but it will not resolve the physics problems. The question of whether the 
QGP is responsible for anomalous J /if) suppression observed by the NA50 
Collaboration will probably and hopefully be resolved in one of the following 
ways (or by combination of some of them): 
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Accurate determination of a a by data on J/ip suppression at E^ = 
160GeV per nucleon. Using this value of a a together with data on 
S(Et) in Pb-Pb interactions might reveal a threshold of "missing con- 
tribution", when plotted as in Fig. || above. 

Combining data on anomalous J/ip suppression with information on 
other signatures, including the data |23| on transverse momenta of sur- 
viving J ftp's and their analyses |2<|, |25|, |2^] and possible future data 



on nucleon number dependence on anomalous J/ip suppression. In this 
respect see e.g. Ref . [pOfl . The explanation of data on anomalous J/ip 
suppression by QGP will only be generally accepted provided that the 
explanation will lead to experimentally verified predictions concerning 
the onset as function of A,B and Et- 

Additional information on pre-resonance interaction with nucleons, see 
e.g. Refs.0, [L7|, [27], ("inverse experiment" and antiproton interac- 
tions in nuclei). 



• Making use of recent results of lattice calculations [p9|1 , which are be- 
coming more and more accurate and lead to the critical temperature 
for the phase transition of about 173MeV. This represents a constraint 
for models with J/ip suppression by a gas of hadrons, since the energy 
density of the hadron gas higher than €qgp{T c ) indicates the transition 
of the hadron gas into the QGP. 
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